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A general description of weak and strong tunneling fixed points is developed in the chiral-Luttinger- 
liquid model of quantum Hall edge states. Tunneling fixed points are a subset of "termination" fixed 
points, which describe boundary conditions on a multicomponent edge. The requirement of unitary 
time evolution at the boundary gives a nontrivial consistency condition for possible low-energy 
boundary conditions. The effect of interactions and random hopping on fixed points is studied 
through a perturbative RG approach which generalizes the Giamarchi-Schulz RG for disordered 
Luttinger liquids to broken left-right symmetry and multiple modes. The allowed termination points 
of a multicomponent edge are classified by a B-matrix with rational matrix elements. We apply 
our approach to a number of examples, such as tunneling between a quantum Hall edge and a 
superconductor and tunneling between two quantum Hall edges in the presence of interactions. 
Interactions are shown to induce a continuous renormalization of effective tunneling charge for the 
integrable case of tunneling between two Laughlin states. The correlation functions of electronlikc 
operators across a junction are found from the B matrix using a simple image-charge description, 
along with the induced lattice of boundary operators. Many of the results obtained are also relevant 
to ordinary Luttinger liquids. 

PACS numbers: 72.f0.-d 73.20.Dx 



I. INTRODUCTION 

Edge states of quantum Hall liquids have attracted 
continuous attention since their importance was pointed 
out by Halperin.El The only gapless excitations of a two- 
dimensional electron gas on a quantum Hall plateau are 
the edge excitations, since the bulk is an incompressible 
quantum liquid.0 A quantum Hall (QH) state with a sin- 
gle condensate, such as v = 1 or v = 1/3, has a single 
bosonic mode of edge excitations, which can be thought 
of as hydrodynamic perturbations of the Hall droplet Bel 
For these states all excitations propagate in a single di- 
rection (the edge is 'chiral'), the direction of the classi- 
cal E x B drift at the edge due to the confining field E 
producing the edge. The chiral-Luttinger-liquid (xLL) 
model of edge structure predicts that more complicated 
states, such as the integer QHE states with v > 2, have 
multiple modes of edge excitations with generally differ- 
ent velocities of propagation, and that some FQHE states 
have modes traveling in both directions along the edge, 
as verified for the v = 2/3 edge by numerical calcula- 
tionsFIl 

This article studies the possible ways of tunneling elec- 
trons or quasiparticles between one FQHE edge and an- 
other or between an FQHE edge and a metallic or su- 
perconducting contact. Such tunneling is of interest be- 
cause, at weak coupling (junction conductance much less 
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than tunneling experiments provide the most sensi- 
tive probe available of edge properties; at strong coupling 
(junction conductance of order ^-) tunneling between 
simple edges is one of the few examples of a solvable, 
experimentally accessible nonequilibrium interacting sys- 



temJS In the remainder of the introduction we outline 
our results and then briefly review the comparison of the 
chiral-Luttinger-liquid description to existing tunneling 
experiments. 

The first part of this article sets up a framework to de- 
scribe possible tunneling fixed points in chiral Luttinger 
liquids and applies it to a number of examples. The best- 
known example of this type of tunneling is between two 
one-component edges such as v = 1/3. Quasiparticle tun- 
neling at finite temperature across a slight constriction 
in a single v = 1/3 quantum Hall bar becomes stronger 
and stronger as the temperature is lowered (quasiparti- 
cle tunneling is "relevant" in renormalization-group lan- 
guage) until the constriction becomes large and the sys- 
tem can be described as weak electron tunneling between 
two separated v — 1/3 edges. The crossover between 
these two fixed points can also be driven by applying a 
voltage across the junction at zero temperature. The fas- 
cinating, experimentally accessible physics of tunneling 
between quantum Hall edges motivates a general study 
of what fixed points are possible in this system. 

Our approach starts by considering the different ways 
in which two edges can be joined, or one edge can be 
terminated. The class of bosonic boundary conditions 
we consider is not formally complete (the set of all al- 
lowed conformal boundary conditions for even a single 
boson is unknown) but includes all previously treated 
cases plus many others. The requirement of unitary time 
evolution gives a physical restriction on possible tunnel- 
ing junctions or edge terminations. As an example of 
the results, it is possible for av = l quantum Hall edge 
to join smoothly to a v = 1/9 edge in the presence of 
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a superconductor (acting as a charge reservoir), but not 
to a v = 1/3 edge: joining a v = 1 and ^ = 1/3 edge 
requires an entropy flow to the surroundings. The ex- 
perimentally relevant correlation functions for electron- 
like operators across the junction can be calculated from 
this framework, and have a simple "image-charge" de- 
scription when the tunneling problem is folded onto the 
half-line. We present a solvable model in which inter- 
actions between edges give a continuous variation with 
interaction strength of effective tunneling charge and the 
I—V curve. This suggests that once interactions between 
edges are considered, nonuniversal features can appear in 
tunneling properties. 

The last section develops a perturbative RG analy- 
sis of the effect on the edge interaction matrix V, de- 
fined below, of random hopping of quasiparticles between 
edge modes. The tunneling properties of an edge are af- 
fected by V for both weak and strong tunneling, and 
for weak tunneling such random hopping (introduced to 
model impuritjLjscattering) is required to obtain univer- 
sal behavioroO. The agreement discussed below be- 
tween the chiral-Luttinger-liquid and composite-fermion 
approaches requires such hopping to drive the edge from a 
nonuniversal starting point to a universal strong-coupling 
fixed point. The importance of V for tunneling motivates 
a study of what form V should take in real systems, where 
some degree of impurity scattering is always present. 

We find that in all principal hierarchy states, random 
hopping of quasiparticles drives the V matrix to fixed 
points where charge and neutral excitations are decou- 
pled. Without random hopping, in general all eigen- 
modes carry charge. When the neutral modes have the 
same velocity, as in IQHE and main-sequence FQHE 
v = n/ (2n ± 1) edges, each edge has one fixed point and 
these fixed points are exactly selpible, as shown by Kane, 
Fisher, and Polchinski (KFP)J§t2l and have higher sym- 
metry than the generic clean system ( "symmetry restora- 
tion by disorder"). Edges with neutral modes in both di- 
rections such as v — 5/7 can have several charge- neutral 
separated fixed points, some of the solvable type found 
by KFP and some unsolvable at present. The RG flows 
have a simple description in terms of the "boost" and 
rotation coordinates for multicomponent %LLs. The RG 
flow equations can a. 
"chiral metal" phase 
uid problems. 

We review one class of experiments to illustrate the 
current status of different descriptions of tunneling into 
edge states. The tunneling I-V exponent from a Fermi 
liquid to .a» edge state is a sensitive measurement of edge 
structural^ (Fig. [I]) which has received much attention 
since the results of Grayson et alx3 showing deviations 
from the predicted curve for filling fractions 1/3 < v < 1. 
Other measurementSLj do show a broad plateau in rough 
agreement with theory, although the plateau is displaced 
from its expected location, perhaps because the density 
profile at the edge is not simply a sudden decline.E-3 
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FIG. 1. The tunneling exponent / oc V a for quantum Hall 
states in the range 1/2 < v < 1. The solid line is the predic- 
tion of the compressible-state theory of Shytov, Levitov, and 
Halperin 16 for infinitely many channels. The dotted line is 
a = 1/v. Solid circles are the main-sequence edges with one 
phase per edge, and other shapes describe fixed points of var- 
ious symmetry classes in edges with multiple phases. 11 The 
states shown are all principal hierarchy states up to 4th level: 
the main-sequence states plus v = 8/11, 5/7, 12/17, 8/13. 

In general, an edge can have multiple edge phases, 
leading to different tunneling exponents. Thus we can 
use tunneling to study transitions between those edge 
phases. Note that for each edge with multiple phases, ex- 
actly one phase has an electron tunneling exponent lies 
on the curve predicted by the composite-fewnian (CF) 
approach of Shytov, Levitov, and HalperirJl3'EZl in the 
limit of an infinite number of channels. Assuming this 
phase is the stable one at finite temperature, the x^L 
and CF approaches appear to include the same essential 
physics, even though the former assumes an incompress- 
ible bulk state and the latter a compressible bulk state. 
The currently unsolvable fixed points are relevant for ex- 
periments (e.g., the phase on the CF line at v = 5/7) 
and are distinguished by having relevant disorder opera- 
tors coupling neutral modes in opposite directions. This 
suggests that the solvability of some fixed points is a 
consequence of their chiral nature. It is worth point- 
ing out that Fig. ^ shows that smoothly varying physi- 
cal properties (in this case, the tunneling exponent) can 
emerge even though the number of edge modes changes 
discretely. 

The outline of the paper is as follows. Section II sets 
up a framework for describing the fixed points of a gen- 
eral quantum Hall edge with a junction or terminus. The 
description can be carried out without reference to the 
details of the microscopic action except for two essential 
formulas, which we derive in section V. Section III applies 
the framework to a number of cases of interest, obtain- 
ing the critical points for a quantum Hall edge coupled to 
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a superconductor and two coupled quantum Hall edges. 
Section IV studies the effects of interactions on tunneling 
in a solvable model, where interactions lead to a contin- 
uous change in the "effective" (i.e., observed in tunnel- 
ing experiments) charge and filling fraction. Section V 
derives the correlation functions and selection rules for 
boundary operators in a chiral Luttinger liquid, which 
can be understood from an "image-charge" picture as in 
electrostatics. Section VI examines the effect of random 
hopping on the interaction matrix of a chiral Luttinger 
liquid, including charge-neutral separation at the edge, 
equilibration of velocities, and the basin of attraction of 
the fixed points of Fig. [I]. Section VII contains a brief 
summary of our main conclusions. 

II. GENERAL DESCRIPTION OF CRITICAL 
POINTS 

A. Characterization of an Abelian edge 

In the following we will give a general description of 
edges of Abelian quantum Hall states, emphasizing the 
data needed to characterize such edges. An Abelian edge 
state with k branches is described by k bosonic fields </> , 
a = 1, .., k. The edge is characterized by a fc-dimensional 
symmetric real matrix K, a fc-component charge vector 
q, a lattice r c that determines the allowed charge excita- 
tions, and a fc-dimensional symmetric real matrix A that 
describes the interaction between edge branches. The 
total electric charge density p e and current j e can be ex- 
pressed through (f> a and q: 

Pe = qaPa, Pa = ^—d x (j) a 

Ztt 

je = qaja, ja = -^dt^a- M 

The allowed charge excitations are created by the vertex 
operators 

v n = e m -*, ner, (2) 

Some V n add charge to the edge and others transfer elec- 
trons or quasiparticles between different edge branches 
but are neutral overall. The charge and the statistics of 
V n are given by 

Q = n T if _1 q, 9 = 7rn T A'" 1 n. (3) 

The vectors in r c satisfy 

i^K^ 1 ^ = integer, 
n,n'er c . (4) 

Thus the charge excitations in T c are bosons or fermions 
with trivial mutual statistics. These vertex operators can 
appear in the Hamiltonian if they have bosonic statistics. 



In many cases, the lattice T c is generated by electron op- 
erators, so we will refer to r c as the E-lattice (although 
r c sometimes contains operators that transfer quasipar- 
ticles between different branches). 

The quasiparticle operators are also labeled by points 
in a lattice, T q : 

T q = {n|m T A" _1 n = integer, for all m 6 T c } (5) 

We will call this lattice the quasiparticle lattice or Q- 
lattice. Note that the Q-lattice T q is the dual lattice of 
the E-lattice T c . From the definition of the E-lattice, we 
see that r c C T q . Since we are going to discuss many 
lattices in this paper, we find it is convenient to use a 
matrix to describe a lattice. We will say a lattice T is 
described by a matrix M if the column vectors of the 
matrix generate the lattice. We will denote such a lattice 
as r = Latt(Af). Also we will use WT to denote the 
transformed lattice of T by W: WT = Latt(WM). Under 
this notation, we can write the E-lattice r c = Latt(C), 
where the k by k matrix C satisfies C T K~ 1 C = integral 
matrix. The Q-lattice can be found to be 

T 9 = Latt{K (C T y 1 ). (6) 

The scaling dimensions of all vertex operators (electron 
or quasiparticle operators) V n can be determined from a 
single matrix A, which depends on K and interaction 
strengths between edge branches: 

h(n) = ln T An. (7) 

Our reason for emphasizing A here rather the "velocity" 
matrix V appearing in the xLL action 

So = J dxdt [K i:j d x 4>idt(t)j + Vijd x <j)id x <j)j] (8) 

is that the universal properties of the junction depend on 
scaling dimensions, determined by A; V contains addi- 
tional information (e.g. the velocities-of the eigenmodes) 
which is unnecessary for this sectionO For example, the 
(equal-space) correlation of V n has a form 

(V n (x,t)Vl(x,0))~r 2h ^ (9) 

The matrix A is not an arbitrary real symmetric ma- 
trix. If we introduce through: 

^ = </ 2 ( / o + \ )Ki/ 2 , (10) 

then A can be expressed as 

A = A7/ 2 5 s - 2 (Af /2 r 1 (11) 

where B st is the boost matrix introduced i n0 and re- 
viewed in section VI: 

B st = exp(® T J) (12) 
and b is a real k + x fc_ matrix. 



3 



B. Termination of an Abelian edge 

First, let us consider the general problem of termina- 
tion of an Abelian edge at a point. We know that a ID 
electron gas contains two branches, one right-moving and 
one left-moving. Such a ID system can be terminated at 
a point (say x = 0) and we can have a system on a half- 
line < x < oo. In contrast, the edge state of a v = l/m 
QH state contains only one right-moving branch. Such 
an edge state cannot be terminated at any point without 
violating unitarity (alternately, the Hamiltonian would 
be non-Hcrmitian) . Now the question is when a generic 
Abelian edge described by (A, q,r c ,A) can be termi- 
nated at a point, and if the edge (A, q, r c , A) can be ter- 
minated, how to characterize the different ways in which 
the edge terminates. Our motivation for considering ter- 
mination of an edge is that a tunneling junction in an 
edge can be considered as a special kind of termination, 
as in the next subsection. 

We find that the edge (A, q, r c , A) can be terminated 
consistently (section V) if there is a 2k by k matrix B 
that satisfies 

B T R- 1 B = 0, 
det(B T B) ^ 

Y B = Latt(B) c fco (13) 

where k = dim(A)/2 and 

f c0 = {n|n e f c , n T K~ l n = even} (14) 

(i.e., the points in T c q describe the bosonic vertex oper- 
ators). The physical meaning of the first two conditions 
is that there are k vectors of length 2k (the columns of 
the matrix B) which are null in the indefinite quadratic 
form A' -1 , orthogonal in A -1 , and linearly independent. 
Previously HaldaneEj described charge-neutral null vec- 
tors of A -1 as "topological instabilities," which allow 
oppositely directed edge modes to localize each other 
and drop out of the low-energy theory; thus condition 
( |l3| ) is that there be k independent topological instabil- 
ities orthogonal in A -1 . (We have relaxed the condi- 
tion of charge-neutrality in order to include situations 
involving coupling to a superconductor.) In order for B 
to exist, A must have the same number of positive and 
negative eigenvalues (the edge has the same number of 
right-moving and left-moving branches). 

If more than one B exists, then the edge (A, q, T c , A) 
can be terminated in more than one way. In other words, 
the boundary at x — can have more than one fixed 
point. These different fixed points will be referred as 
different terminations of the edge. 

For a termination labeled by the matrix B, the fields 
4> satisfy the following boundary conditions at the termi- 
nation point 

B T (j> = 27rn, n = integral vectors. (15) 



Note that <f> can satisfy different boundary conditions (la- 
beled by n) even for a single type of terminate labeled 
byB. 

The allowed charge excitations (vertex operators) at 
the boundary are labeled by points in a fc-dimensional 
lattice T qB (called the boundary Q-lattice): 

r 9B = (16) 

Latt [kB{B T B)- 1 - ^B(B T B)- 1 B T KB{B T B)- 1 ^ 

The boundary quasiparticle operators have a form = 
e dT< ^, 1 € r gB . The scaling dimension of is given by 

h b {l) = \ T k- 1 B{B T KB)- 1 B T K- 1 \ (17) 

This is one of the main results of this paper. 

In the above discussion of termination points, we have 
ignored any symmetry properties and the related selec- 
tion rules. In particular, the boundary condition char- 
acterized by B may not conserve electric charge. As a 
result, a boundary vertex operator may not carry a defi- 
nite electric charge. In order for the termination labeled 
by B to conserve the electric charge, we must require the 
B matrix to satisfy 

B T K~ 1 q = (18) 

For the charge conserving termination points, the electric 
charge of a boundary vertex operator is found to be 

Q = q T k-H (19) 

For a general termination described by B, there are k 
combined charges that are conserved near the boundary. 
Their densities are given by 

Pa = B ba d x 4> a /2rr (20) 

The boundary operator carries definite values of these 
k combined charges: 

Q = B T A- 1 1 (21) 



C. Tunneling junction in an Abelian edge 

Now let us consider a more practical problem - a tun- 
neling junction in an Abelian edge. Assume the edge is 
described by (A, q, r c , A). The tunneling is described by 
vertex operators in the Hamiltonian, V n , n £ T c (the lat- 
tice F e determines which vertex operators are allowed). If 
the charge is conserved, only the neutral vertex operators 
in r c can be added to the Hamiltonian to describe tunnel- 
ing. These neutral vertex operators describe the different 
types of electron/quasiparticle tunneling between differ- 
ent edge branches. We would like to remark that the 
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tunneling junction referred to here can also be viewed as 
a defect or an impurity on the edge. 

To describe the possible quantum fixed points of the 
tunneling junction, we can fold the edge in (—00, 0) on 
top of the edge in (0, 00) by introducing 2k fields <p a on 
(0,oo): 



4> a {x) = (f>a{x) 
4>k+a(x) = ~(j)a{-x) 

x > a = 1, .., k. 
The resulting edge is described by (K, q, r c , A): 



(22) 



K = 
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A = 
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q = 
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!■ 



(23) 



The edge is terminated at x — 0. Now the problem of the 
different fixed points of a tunneling junction becomes a 
problem of different ways that the edge (K, q, r c , A) can 
terminate at x = 0. More precisely, each fixed point of 
the tunneling junction correspond to a way in which the 
edge (K, q, r c , A) terminates. 



D. Joining of two Abelian edges 

Next, let us consider the problem of possible different 
edge states for a given QH liquid. There is the possibil- 
ity that different edge potentials and/or electron interac- 
tions can lead to different edge states without changing 
the bulk QH liquid. Indeed, it has been shown that at 
zero temperature there are multiple stable edge phases 
of some disordered FQH states with neutral modes in 
both directionsliil, such as v = 5/7. If a given QH liquid 
does have different edge states, then we can put different 
edge potentials on different segments of the edge, leading 
to different edge states on different segments. Thus two 
different edge states of a given QH liquid can always be 
connected together. This motivates us to ask the follow- 
ing question: when can we connect two Abelian edges 
(Ki,qi,Tcx, Ai) and (K 2 , q2, r c2 , A 2 ) at a point x = 0? 

Again we can fold the edge (K 2 , q 2 , r c2 , A 2 ) in (—00, 0) 
on top of the edge (K\, qi, F c i, Ai) in (0, 00). The result- 
ing edge is described by (K, q, r c , A): 



K 



A = 



Ki 

-K 2 

Ai 

A 2 



r c = r cl e r c2 

qi 
q2 



The two edges (K\, qi, T c i, A x ) and (K 2 , q 2 , T c2 , A 2 ) can 
be joined together only if the edge (X,q, T c , A) can be 
terminated. The different ways to join the two edges 
correspond to the different ways to terminate the edge 
(K,q,f c ,A). 



III. SOME SIMPLE EXAMPLES 

To gain a more intuitive understanding of the results 
summarized above, we would like to discuss a few simple 
examples. 



A. v — 1/m edge state coupled to superconductor 

First let us consider an edge of a v = 1/m Laughlin 
state. We place a tunneling junction to a superconduct- 
ing state at x = 0. After the folding, we get a two-branch 
edge: 



/ m 









\ 


— m 1 


, q = 






K 



_ f 1/m 
^ ~ 1 1/m 



(24) 



To understand 



The E-lattice T c is generated by n T = (m, 0) (which 
creates an electron) and n T = (1,1) (which transfers a 
quasiparticle of charge 1/m from one edge to the other). 

That is, f c = Latt(C), C = (™ J 

the fixed points of the tunneling junction, we first study 
the termination of the two-branch edge. 

One termination (called fixed point A) is described by 

B = ( J ^ . Such a termination conserves electric charge, 

since B satisfies (|l8|). To obtain the boundary quasipar- 
ticle operators, we need to find the boundary Q-latticc 
T qB , which is given by (p"7|). We find 



An element 1 of the lattice T qB has the form 



1 = I 



m/2 
-m/2 



(25) 



(26) 



where / € Z . Thus the boundary quasiparticle operator is 
labeled by a single integer We denote such an operator 
by Vf. The scaling dimension of V l is 

h\l) = l T K- 1 B(B T AB)- 1 B T K- 1 l = ml 2 /2, (27) 

and the charge is Q = c^K^ 1 ! = I. It is clear that the 
boundary vertex operator V± just creates an electron. 
It is also clear that Vi cannot appear in the boundary 
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hamiltonian, since only pairs of electrons can be added 
or subtracted from the superconductor. 

To determine the boundary operators that can appear 
in the Hamiltonian, we need to consider charge conserva- 
tion. In this problem the charge (Q mod 2) is conserved. 
Thus only Vf with even I can appear in the boundary 
Hamiltonian. The leading operator that can appear in 
the boundary Hamiltonian is V% j which has scaling di- 
mension h = 2m and a charge Q — 2. V% adds two 
electrons to the boundary and describes the tunneling to 
the superconductor. The scaling dimension of V% deter- 
mines the stability of the fixed point A. Since the scaling 
dimension 2m > 1, the fixed point is stable. Physically, 
the fixed point A corresponds to a junction in which the 
tunneling to the superconductor vanishes at low energies. 

Another termination (called fixed point B) is described 



by B = 



Such a fixed point does not conserve 



the charge. The boundary Q-lattice T qB is 



The elements 1 in r o 



Latt( ^ -|y 2 
have the form 



\-l( 1 ' 2 
1 1/2 



(28) 



(29) 



and the boundary quasiparticle operator is labeled by an 
integer I. The scaling dimension of V/ 3 is 



h b = I 2 /2m. 



(30) 



Note that although the charge Q — Qi + Q2 is not con- 
served, the charge difference Qd = Qi — Qi is conserved 
at low energies. The Qd charge of VJ b is given by 



l/m 



(31) 



where qj = (1,-1). This means that the boundary oper- 
ator Vi transfers l/m charges between the two branches. 

The leading operator that can appear in the boundary 
Hamiltonian is V b , since we allow any amount of charge 
to be transferred between the two branches. Its scaling 
dimension is h b = l/2m. Since l/2m < 1 the fixed point 
B is unstable. Physically the fixed point B corresponds to 
a junction with strong tunneling to the superconductor. 
A low energy incoming electron will be scattered into an 
outgoing hole by the junction. 

Let us start with a v = l/m edge state coupled 
strongly to a superconductor at a point. At high energy 
scales (temperature or applied voltage), the junction is 
close to the fixed point B. As the energy is lowered, the 
coupling strength flows to zero and the junction flows 
from the unstable fixed point B to the stable fixed point 
A. Note that the scaling dimensions of leading boundary 
operators at the fixed points A and B are given by 2m 
and l/2m, which are inverse of each other. Thus the 
fixed point A and B form a duality pair. The crossover 



between fixed points A and B can be solved exactly and 
has been studied extensively in tunneling between FQH 
edge statesB. In section IV we show how a simple model 
incorporating interactions between electrons on different 
edges can be mapped onto this exact solution. 



B. (y\,v-x) = (1/mi, 1/7712) edge state coupled to 
superconductor 

Second, let us consider a two-branch edge described by 

A superconductor covers and couples to the edge in the 
region (— 00, 0). We would like to ask: can the coupling 
produce a gap in the (— 00, 0) region, or equivalcntly: 
can the edge terminate at x = 0. We find that if ^ is 
a square of a rational number, then the B matrices exist 
and the edge can be terminated. These edges include 
(mi,m 2 ) = (m,m), (1,9), (3,27), (9,25), etc. Note that 
the (mi, ma) = (1,9) edge is the edge of the v = 8/9 QH 
state. Thus the gapless edge excitations of a v = 8/9 QH 
state can disappear when coupled to a superconductor. 

First, we discuss the case (7711,7712) = (m,m). This 
case was discussed above, except that the lattice T c is 
different for the present setup. Even here there are two 
different possibilities. In the first setup where the two 
edges belong to the same QH liquid (Fig. |a), f c is given 
by 



f c = Latt(^ 1 



1 m 




), 



(33) 



since we can transfer quasiparticles of charge l/m be- 
tween the two branches of the edge. This setup is identi- 
cal to the case discussed in the last subsection, and here 
we just repeat the results obtained before. The fixed 

r 



point A is described by B = 



The leading bound- 



B = 



The leading boundary operator that can 



ary operator that can appear in the hamiltonian has scal- 
ing dimension 2m. The fixed point B is described by 
m 
— m 

appear in the hamiltonian has scaling dimension l/2m. 

In the second setup where the two edges belong to two 
different QH liquids (Fig. |b), r c is given by 



r c = Latt((™ -]) 



(34) 



since we can only transfer electrons between the two 
edge branches, not quasiparticles. Here we can define 
two charges. The first one is the total electric charge Q 
described by the charge vector q. Only (Q mod 2) is 
conserved. The second one is the difference of the elec- 
tric charge on the two branches Qd = Qi — Qi- Qd is 
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described by a second charge vector q,j 
(Qd mod 2) is conserved here. 



Only 



(a) 








(b) 




FIG. 2. Two schematic setups for a superconducting con- 
tact at the termination of two edge states. In (a) the edge 
states belong to the same droplet, in (b) to separate droplets. 



charge Qd mod 2 = mod 2, and has scaling dimension 
2m. 

In the above, we see that not every boundary operator 
in r& can appear in the Hamiltonian. One may wonder 
what is the meaning of the operators which cannot ap- 
pear in the Hamiltonian. We notice that although there 
are no gapless excitations in the region (— oo,0), there 
are degenerate ground states in that region. Thus there 
are domain- wall-like excitations in (— oo, 0) region. The 
ground state changes from one to another as we go across 
a domain wall. A domain-wall-like excitations is labeled 
by a vector 1 in the boundary Q-lattice T QB . It carries 
the same quantum numbers as the boundary vertex op- 
erator labeled by the same vector 1. As we move a 
domain-wall like excitation to the boundary, it becomes 
a boundary quasiparticle excitation described by Vf. 



f 



(a) 



v = I/ni 



1 SC , 



The fixed point A is described by B 
boundary Q-lattice r„ B is 



T qB — Latt( (_y2 



The 



(35) 



The boundary quasiparticle operator is labeled by / with 



1 = / 



1/2 
-1/2 



. The scaling dimension of V/ 3 is 



h b = I 2 /2m, 



(36) 



and total charge Q is conserved at low energies. The 
boundary quasiparticle operator carries definite 
charge Q = K _1 \ — l/m. The leading boundary 
quasiparticle operator that can appear in the boundary 
hamiltonian is V^m which carries charge Q mod 2 = 
mod 2, and has scaling dimension 2m. 

The fixed point B is described by B = ( ) . The 



boundary Q-lattice T qB is 

r te =Latt(Q/^) 



(37) 



The boundary quasiparticle operator is labeled by I with 
1 = Z [ 1/2 ) - ^ e scann § dimension of VJ b is 



h b = I 2 /2m, 



(38) 



For the fixed point B, the charge Qd is conserved at low 
energies. The boundary quasiparticle operator V/ 3 carries 
definite Qd charge: Qd — c^K~ x \ = l/m. Such an op- 
erator transfers I /2m charges between the two branches. 
The leading boundary quasiparticle operator that can ap- 
pear in the boundary hamiltonian is V 2 b TO which carries 



(b) 



v-l.'ni 



v-l.'ni 



FIG. 3. Two schematic setups for a superconducting con- 
tact to two edge states. In (a) the edge states belong to the 
same droplet, in (b) to separate droplets. 

To physically measure the effect of the boundary quasi- 
particle excitation, let us consider the two setups in Fig. 
H where only a segment of edge is covered by supercon- 
ductor (The two setups in Fig. || are related to the two 
setups discussed above). In those setups, we can create a 
pair of boundary excitations: the first one is created by 
Vj ', 1 <E T;, at one boundary and the second one is created 
by V_i at the other boundary. Such an operation trans- 
fers a domain-wall-like excitation from one boundary to 
the other, and is an allowed process. For the fixed point 
A in the first setup, the lowest scaling dimension for the 
above pair-creation operators is m. For the fixed point 
B in the first setup, the lowest scaling dimension for the 
pair-creation operators is l/m. For m > 1, the fixed 
point A is stable while the fixed point B is unstable. For 
both fixed point A and B in the second setup, the lowest 
scaling dimension for the pair-creation operators is l/m. 
For m > 1, the operator is relevant (since l/m < 1), and 
the two edge states separated by the superconductor will 
join together at low energies. 

Next, we would like to discuss the (mi, 7712) = (1,9) 
case. We will also consider the possibility of short-ranged 
density-density interactions between modes on the mi = 
1 and on the 777,2 = 2 edge. The interaction effects being 
discussed are described in detail in the next section for 
a different case; to avoid duplication we will simply start 
here from the boost form ( |ll|) of the matrix A: 
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1 
1/3 



/ coshr sinhr 
I sinh r cosh r 



1 

1/3 



(39) 



The boost parameter r measures the interaction 
strength: r = corresponds to unmixed v = 1 and 
v = 1/9 states. The lattice for charge excitations T c 
is given by 



f c = Latt( 



1 
9 



(40) 



The first termination (called fixed point A) is described 
9 



by B = 



The boundary Q-lattice is given by (|l7|): 
= Latt ((-l/ 2 )) (41) 



r 



The boundary quasiparticle operators are labeled by an 
integer / with 1 = I ' ^' ^ 



The boundary operator 



-1/2 

V b has scaling dimension h b 

combined charge <5(i.3) = Qi + described by the 
vector ^ = (1,3) is conserved at low energies. The 
Q(i,3) charge of V b is Q(i,3) = cff^K^l = 1/3. The 
leading operator that can appear in the Hamiltonian is 
Vq which has an even number of <5(i,3) charges and a 
scaling dimension h b = coshr — sinhr. 

The second termination (called fixed point B) is de- 
scribed by a lattice B = ( ^ ] . The boundary Q-latticc 



|g (cosh r — sinh r) . The 



is given by 



T qB = Latt( 



1/6 
1/2 



)• 



(42) 



The boundary quasiparticle is labeled by an integer I 
1/2, 



with 1 



The boundary operator V b has 



scaling dimension h b = (coshr + sinhr). The com- 
bined charge (5(1,-3) — Qi — 3Q2 described by the vec- 
£ _ 31 = (1,-3) is conserved at low energies. The 



tor q£ , 



<5(i -3) charge of V; is Q(i _ 3 ) = q£ _ 3 )i£' T = £/3- The 
leading operator that can appear in the Hamiltonian is 
Vq which has an even number of Q(i _ 3 ) charges and a 
scaling dimension h b = coshr + sinhr. Note that since 
(coshr — sinhr) (coshr + sinhr) = 1, fixed points A and 
B form a duality pair. 



C. Junction between v\ = l/mi, ^2 = l/wi2 edge states 

This case has been treated previously in-the. absence of 
intermode interactions by several authorscJoEicj. Here 
we will show how the known results for the two fixed 
points are recovered in our framework, and in the next 



section consider interaction effects. A new result even 
in the absence of interactions is the information found 
below about the lattice of charged boundary operators. 
For definiteness we consider the case of tunneling between 
v = 1 and v = 1/3 states which has attracted the most 
attention. At weak coupling between the two edges, the 
most relevant neutral operator tunnels an electron be- 
tween the two edges, with scaling dimension 2; at strong 
coupling the most relevant neutral operator has scaling 
dimension i. This neutral operator can be interpreted 
as tunneling between different minima of the boundary 
cosine interaction in the sine-Gordon model, or as quasi- 
particle tunneling in an effective model of two v = \ 
edges. 

With no intermode interactions, we have 



K = 












)... 










(43) 



Fixed point A (weak coupling) is described by 



B 




(44) 



The meaning of this B matrix is that the incoming and 
outgoing v = 1 edges are joined continuously to each 
other, and similarly for the incoming and outgoing v = 
1/3 edges. The boundary Q-lattice is 



— Latt( 



( I 


\ 


1 

2 








1 

2 


V 


-V 



(45) 



A general boundary operator 0\ with 1 = 4,0,0) 
+^(0, 0, |,— |), l\ and I2 integers, has electric charge 
Qi = h + h/2 and scaling dimension Zi 2 /2 + ?2 2 /6. The 
most relevant neutral boundary operator has (£1,^2) = 
(1,-3) and scaling dimension 2, as expected for electron 
tunneling at weak coupling. 

Fixed point B (strong coupling) is described by 



B 




(46) 



The meaning of the first column of B is that two incoming 
electrons become one outgoing electron and three outgo- 
ing q = i quasiparticles. In the boundary sine-Gordon 
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language, this corresponds to pinned (Dirichlet) bound- 
ary condition on the neutral mode. The second column 
is simply the charge vector, indicating that overall charge 
is conserved across the junction. The boundary Q-latticc 
is 



Lattf 



/ 



_3_ 

10 

10 
7_ 

10 
_7_ 

10 



10 

_4 
5 

3 
2 
.! 



(47) 



A general boundary operator 0\ with 1 = 
l i(jo> To>-To>~To) +l 2(-To>~i> I' f) has charge Q = l 2 
and scaling dimension h 2 /2 — 3M2/2 + 3?2 2 /2. The most 
relevant neutral boundary operator has — (1,0) 

and scaling dimension ^ , as expected. 



IV. EFFECTS OF SHORT-RANGE 
INTERACTIONS 

In this section we consider the effect of short-range 
interactions on tunneling through a point contact be- 
tween two Laughlin states. In the absence of interac- 
tions, the nonlinear I — V c.urve was found exactly by 
Fendley, Ludwig, and SaleurO via a mapping onto the in- 
tegrable boundary sine-Gordon (BSG) model. First we 
show in the BSG formalism that a simple solvable model 
incorporating interactions gives a continuous renormal- 
ization of the effective fractional charge appearing in the 
I — V characteristic. We use the BSG formalism since we 
will eventually be interested not only in the fixed points, 
which can equally well be described in the B-matrix for- 
malism of section III, but also in the crossover. The I — V 
curve measured in tunneling experiments in real systems, 
where the screened Coulomb interaction is present, will 
thus be sensitive in some geometries (discussed below) to 
nonuniversal electron-electron interactions. Our model 
is different from that of Pryadko et alE3, which uses a 
long-range Coulomb interaction regularized by an open- 
ing angle at the junction. 

The effective action describing tunneling between 
edges of two Laughlin states with filling fractions v\ = 
1/mi, V2 = l/ni2 is 

S = Sfroc + >5tuni 

SWee = J dxdt y^l\,.<l,o,<),Q. - Vijd x cf>id x (f)j], 



S tua = T5{x){e imi ' l > 1 - im2<f > 2 ). 



(48) 



Here the matrix K, which describes the statistics of the 
vertex operators e lUi ^ i created from the bosonic fields, is 



K =(7/ ° ) 

— 7712 



(49) 



where we have taken the two edges to propagate in oppo- 
site directions. If V is diagonal the physics is independent 



of whether the modes are copropagating or counterprop- 
agating, but we are interested in the case of general V 
in which case there are differences, as seen below. If the 
positive definite matrix V is diagonal, its two entries V\\ 
and V22 are the velocities of the two modes. The off- 
diagonal elements of V correspond to a density-density 
interaction across the two edges, since the electron den- 
sity is proportional to d x <pi for each mode. 

The above action maps onto a boundary sine-Gordon 
model, with boson radius determined by the filling frac- 
tions of the original states and by the matrix V. The 
boundary sine-Gordon model contains one nonchiral bo- 
son (i.e., with both left and right components) on the 
half-line. The mapping consists of rotating the fields 
4>i , 4>2 so that one new combination <f>\ is proportional 
to the exponent m\(j>i — m2<p2 in <5tun, while <f>2 does not 
appear in Stun and hence is free. Then folding the field 
<t>\ onto the half-line and rescaling gives the action 

(50) 

The constant f3, given for diagonal V by (5 = 

^1 + ^L), gives the tunneling term in (|5(]) the same 

scaling dimension A = (mi 2 + TO2 2 )/2 as in (|48"|). The 
velocities of the edge modes, defined as the velocities in 
a basis where V is diagonal, should strictly speaking be 
equal for this rotation of fields to be valid, but since the 
tunneling takes place at a point and there is no coherence 
along the edge, a difference in velocities should not have 
much effect. 

In order to calculate the conductance across the tun- 
neling junction, the effective f3 which appears in 5bsg 
needs to be determined, as well as the contribution q e g 
to the current from each tunneling event. Previously only 
certain discrete values of /3, corresponding to tunneling 
between Laughlin states, were thought to be physically 
relevant for edge tunneling. This is because a general (3 
describes tunneling between two chiral Luttinger liquids 
with continuous Luttinger parameter, but only specific 
values of the Luttinger parameter correspond to quan- 
tum Hall states v = 1/m. The main result of this sec- 
tion is that tunneling between Laughlin states with non- 
diagonal V is described by the boundary sine-Gordon 
model with continuously varying [3 and q c R . 

The model which we solve exactly has a region of con- 
stant interaction strength (between contacts V\ and V 2 in 
Fig. |j) and zero interaction elsewhere. It is essential that 
the two modes in the interaction region be oppositely di- 
rected, so that the scaling dimension of the tunneling 
operator is affected by V. The first step is to write the 
positive definite matrices V and A in terms of a "boost" 
parameter r. The advantage of doing so is that A is only 
a function of r and not of the eigenmode, velocities Vi 
which affect V; the boost decompositionlid isolates the 
dependence of A on as few parameters as possible. 







V = 


K^ 2 B 


A = 


K^B 


K^ 2 = 


( y/rni 

\ o 


B = 


/ cosh- 


\ sinhi 



v 2 



BK 



1/2 



(51) 



Now the scaling dimension and transferred charge of the 
tunneling operator can be simply expressed in terms of 
A(r). In the following we will specialize to the case 
mi = 1, rn2 = 3, i.e., a v\ = 1 edge and v 2 = \ edge 
propagating in opposite directions. The results generalize 
simply to other values of m\ and m 2 , although the bound- 
ary sine-Gordon crossover result used below requires that 
at most one relevant operator be present, restricting m\ 
and 77i2 somewhata. 



Vi 



Vl 



V 2 



FIG. 4. Possible experimental geometry for point tunneling 
between quantum Hall states v\ and V2- The density-density 
interaction between edges is nonzero in the shaded region, and 
zero elsewhere. The two contacts at voltages Vi and V2 are 
assumed to populate edge modes propagating away from the 
contact up to energy eV. 

The intuitive meaning of the effective charge transfer 
can be understood by considering the "charge-unmixed" 
point t = where one of the two eigenmodes (i.e., modes 
which diagonalize K and V) is neutral, and one charged. 
Then in this basis the tunneling operator, which is neu- 
tral, is of the form exp(iC<f> n ) + h.c. for some constant 
C . The tunneling operator for this value of r does not 
transfer any charge from one eigenmode to the other, 
since only one eigenmode carries charge. Hence the con- 
ductance measured across the junction is independent of 
the rate of tunneling events determined by the coefficient 
of the tunneling operator. The value r = log(l + \/3/\/2) 
corresponds to decoupled v = 1 and v = 1/3 edges in the 
interaction region, idiida-|has been previously studied in 
a number of works .II3EIE3 

The scaling dimension of the tunneling operator 
exp(zmj</>i) is imA(r)m. The effective tunneling charge 



is determined by how far the neutral tunneling opera- 
tor m is from being an eigenmode of the system: q e g = 
|qAm, which is zero if m is an eigenmode and grows as 
A moves away from the charge- unmixed point. The cur- 
rent measured across the two contacts of Fig. |] if there 
is no electron tunneling across the junction (if V\ — V> is 
small) is I — (Vi — ^)(| + + 02)), where o\ and 02 
are the conductances along the edges in the interaction 
region of Fig. ^. We use the identity <j\ + 02 = tAt in 
what follows. 

The total change in conductance from small V to large 
V is fixed by the scaling dimension of the electron tun- 
neling operatcp,and the effective tunneling charge. From 

previous worked, it is known that if the effective tunnel- 

2 

ing charge is 1, the conductance change is j^, where A 
is the scaling dimension of the electron tunneling oper- 
ator. For instance, in tunneling between two v — 1/3 
states, the conductance change between no tunneling 

(two Hall droplets) and no backscattering (one Hall bar) 

22 

is |r = In the case of interactions, e must be replaced 
by the effective charge transfer per tunneling event q e g. 
Hence for the system of Fig. [|, 



^max 2 

~ 3 
2 



tAt 2 + cosh(2r) 



e 



2 
tAt 



3 

(tAm) 2 



2 2m Am 

sech(2r) 



(52) 



In fact the whole conductance curve between these two 
values can be calculated from the mapping to Sbsg- Be- 
fore doing so, there is a simple check on our results which 
gives some insight into why the above values are natu- 
ral. Assuming conservation of energy (i.e., no dissipation 
at the junctions) gives two possible values of the cur- 
rent: / = (Vi — V2X! + ^(cti + 02)), corresponding to 
no tunneling current, and / = (Vi — VOd + g7^+^j )- 



The two corresponding values of the conductance are 
exactly those in ([32]). Thus our calculation reproduces 
the asymptotic values of the conductance consistent with 
zero-dissipation fixed points. 

The current-voltage characteristic can be calculated 
(Fig. H up to one overall constant in the energy scale, 
which corresponds to the initial strength of tunneling.El 
The result is, with V = Vi — V2 and t7 ma x as in (|5T 



/ = 


CmaxV 




//(I) 


/« = 


9 2 eS V , 



if V > TsA-^-^VA^I ' 



r(2) _ 



hA 



1 - 



^ /n(l/A) \ 



n=l 

(-l)" +1 0F 7 (ng) (eV_ 
2T(n)r(3/2 + n( ff -l)) \T B ' 



2n(g-l) 



(53) 
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Here Tg is some cutoff-dependent constant which may 
vary with the boost parameter r. This calculated current 
should be relevant as long as the interaction strength is 
nearly constant in the region around the tunneling junc- 
tion. The details of the interaction far away from the 
tunneling junction should not matter as long as current 
is conserved in the incoming and outgoing edge branches. 




. 8 



1 2 3 4 5 

C(t)V 

FIG. 5. Differential conductance dl/dV in units of e 2 /h 
versus scaled voltage C(t)V for three different values of in- 
teraction strength (boost parameter) r. The horizontal axis 
is expected to scale by a different cutoff-dependent constant 
C(t) for each value of r. 

V. CORRELATION FUNCTIONS IN BOUNDED 
LUTTINGER LIQUIDS 

A. Review of Abelian edge states 

First we review _the xLL description of the edge of an 
Abelian QH stateia The collective modes that propagate 
along the edge can be described by several density oper- 
ators p a (x), which satisfy the algebra 

[p a (x), p b (y)\ = ^-(K-Xrf'ix - y) (54) 

where K is a symmetric integral matrix that character- 
izes an Abelian edge. The total electric density is a sum 
of p a weighted by a charge vector q a : 

Pe = q a Pa- (55) 

If there is no binding between electrons, we can choose 
the symmetric basis in which q a = 1 . In the following we 
will work in this symmetric basis. 

The edge dynamics is described by the Hamiltonian 

H = J dxTTV a bPa(x)pb(x) (56) 

at low energies, where V is positive definite symmetric 
matrix. The corresponding Lagrangian density is given 
by 



C = -—(K ab <j )a dtd x <p b -V ab <j> a dl<p b ) (57) 

where d x 4> a /2n = p a . 

For a circular edge (of length L), if we assume <j) a to 
satisfy the periodic boundary condition (f> a (0) = 4>a{L), 
then the </> only describe the neutral fluctuations. To 
go beyond the neutral fluctuations, let us first consider 
the charge fluctuations created by electron operators ip e . 
There are dim(A") different electron operators 

i> b e {x) = e l < b ^ x] (58) 

where is the b th column of K: = K ab . An multi- 
electron operator has a form e m < s 't>a( x ) w jth n e belonging 
to a lattice r e , which will be called the electron lattice. 
The lattice T e = Latt(A^), since T e is generated by the 
column vectors of K. 

The charge fluctuations created by the electron oper- 
ators can also be described by the <fi fields in the La- 
grangian. But now the <fi fields no longer satisfy the pe- 
riodic boundary condition. Instead they satisfy 

<t>a{L) ^M0) + 2nN a (59) 

where N a are integers. In fact, the above </> describe 
excitations created by Yibi^e)^ fr' om the neutral ground 
state. 

In general, the (multi-)electron operators do not repre- 
sent all possible "charge" excitations. The most general 
charge excitations are created by 

V> c (x) = e mMx} (60) 

where n c are integral vectors belonging to a lattice T c 
(which will be called the E-lattice). The general charge 
excitations contain all excitations created by the electron 
operators, and thus T e C T c . The general charge excita- 
tions can also include excitations that transfer fractional 
charges between different edge branches, but have integer 
overall charge. 

The E-lattice r c must satisfy certain conditions. Since 
ip c (x) all carry integral charges (which can be zero), the 
vectors in T c satisfy 

q T A" _1 n c = integer (61) 

The operators ip c (x) are also mutually local; that is 

(n / c ) T iT- 1 n c = integer (62) 

for any n c and n c in r c . This condition implies that 
r c C Latt(A" _1 ). Just like the electron operators, the 
excitations created by t/j c = e mc '^ are described by cj) 
that satisfy the boundary condition 

(j){L) = <j){Q) + 2nK- 1 n c (63) 

We can also use the above to define the periodicity of the 
<f> field. That is (j) and <fi' are equivalent if 

4> = 4>' + 2itK~ 1 xi c (64) 

for some n c £ f £ . 
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B. Correlations with boundary 



[B^ a (x),B b j4, b (y)]=0, for all/, J. 



(71) 



Now we add a tunnel junction or termination to the 
system. In what follows it is assumed that in the junc- 
tion case the edge has been folded onto the half-line as 
in section II, so that one set (K, q, T c , A) describes both 
incoming and outgoing edges, with a total of 2k bosonic 
fields. It will be shown below that the number of incom- 
ing fields must equal the number of outgoing fields. In 
the presence of tunneling between different edge branches 
at x = 0, the Lagrangian also contains terms of type 



ti5(x) 



h.c. 



(65) 



where tj is the (real) tunneling strength, and the integral 
vectors n 1 are points in the T c lattice which also satisfy 



(n j yK-W = even. 



(66) 



(Thus e m are bosonic operators.) If electric charge is 
conserved at the boundary, n 1 should also satisfy 



c^K-W = 



(67) 



so that e m represents a neutral operator. For the time 
being, we will not impose the above charge-conservation 
condition. 

The Lagrangian density on the half-line is given by 



~ ((K) ab 4> a dtd x 4> b ~ (V) ab 4>adl&) (68) 



C = — 



We also need to specify what are the allowed charge exci- 
tations, in addition to the Lagrangian (|68|). This can be 
achieved by specifying the periodicity conditions of the 
6 fields: 



<j>' = 4> + 27r^ _1 n c 



(69) 



where n c are vectors in the lattice f c . 

To completely define our system on the half-line, we 
need to include boundary conditions at x = 0. First let 
us consider the boundary condition in Hamiltonian lan- 
guage. We will start with the following type of boundary 
conditions, specified by a collection of vectors Bf with 
I=l,...,k: 



B^ a (0) = 0, for all I 



(70) 



Since we are looking for critical points, we would like 
to study boundary conditions that are invariant under 
scaling, suggesting the form (|70|). However, the above 
form is inconsistent with the periodicity conditions on (p. 
It turns out that (|70|), although sufficient to determine 
the correlation functions of vertex operators, must be 
improved (later in this section) in order to determine the 
lattice of allowed boundary quasiparticle operators. 

Since </> a (a;)'s do not commute with each other, the 
above boundary conditions are self-consistent only if 



This implies that the vectors that characterize the bound- 
ary condition must satisfy 



B T R- 1 B = 0, det(B T B) ^ 



(72) 



where B is the 2k by k matrix formed by Bj. Note 
that two different B matrices, B\ and B2, specify the 
same boundary condition in the sense of ( |7fj| ) if they are 
related by 



B 1 = B 2 U, U E L{k) 



(73) 



where U is an k x k invertible real matrix U . In this case 
we say B\ and B2 are "weakly" equivalent 



Bi ~ B 2 



(74) 



In the Lagrangian language, the above boundary con- 
dition corresponds to B T (j>(0) = on the fields <j>. Only 
for certain choices of B can the operator K, — (K) ab dtd x — 
(V) ab d^. be hermitian. We will show that the condition 
on B that makes JC hermitian is nothing but (|7^). To see 
this, we first note that 



dxdt(f>2 fC(f>i — / dxdttfii K,<j)2 + / dt(dt(j> 2 K4>{) x =o 



Thus JC is hermitian if B is such that 

4> a K ab 4>' b - 



(75) 



(76) 



for all <j), 4>' that satisfy B T tp — B J <jf — 0. Let us choose 
a basis that the vectors in the null space of B have the 
form (fi^uu = (0, 0, a, b), i.e., the first k elements are 
zero.. In this basis, by ( f76| ) K has the form 

A" I "V W ! ■ (77) 



K{ 



and B has the form 



B = 



(78) 



where K\ and B\ are invertible. Since K 1 has the form 



K~ 1 = 





(Kir 1 



(79) 



thus the condition (^6|) implies the condition (|7^) . We see 
that the hermiticity of the operator (K) ab d t d x - {V) ab dl 
also requires B to satisfy (|72|). 

Let us summarize our results so far. The critical 
boundary conditions of an Abelian edge described by K 
are characterized by a matrix B that satisfies (^). If 
such a B does not exist, then the edge state described 
by K cannot be terminated at a point. If such a B does 
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exist, then we can consistently impose a boundary condi- 
tion B T <f) = and terminate the edge state at x = (pro- 
vided that B satisfies some other conditions that will be 
discussed later). If more then one inequivalent B exists, 
then the edge can be terminated in more than one way 
(the tunneling junction has more than one fixed point). 

Now the question is when B exists. First, from ( j77| ) 
we see that the signature of K must be zero in order for 
B to exist. Since K is invertible, we can write K as 



K = K 



1/2 



Iky. k 







-Ikxk 



K 1/2 



(80) 



We see that B always exists for the above K. A generic 
B that satisfies (72) can be written as 



B 



1 

-L n 



Kl/2 



(81) 



Thus B exists if and only if K has a vanishing signature 
(there are as many incoming modes as outgoing modes) . 
In this case the different boundary conditions are labeled 
by an element in the 0(k) group. 

Next we would like to calculate the correlation func- 
tions of the fields 4> in the presence of the boundary. First 
we choose a basis in which both K and V are diagonal 
(this can always be doneta): 



K 



I 
-I 



V 









-VL 



V dl 



(82) 



The block vr has positive diagonal elements which are 
the velocities of right movers, while Vl has negative di- 
agonal elements which are the velocities of left movers. 
Thus the velocity for each branch is given by v a — 
(S3Vdiag)aa- In this basis a generic B is equivalent to 
a simple form: 



B 



rpT 



I 



(83) 



Let Go : ab(x,t) be the correlation of <p a (x.t) and 
4>b(fi) in the absence of the boundary: Go^ a b(x,t) — 
(4> a {x, i)0b(O)^} . The correlation satisfies a linear equa- 
tion 

(-k ac d Xl d t + V ac dt) G , cfc (x, t) = 6 ab 6(x)6(t) (84) 



Since K and V are diagonal in the present basis, the 
above can be rewritten as 

^ (-K aa d x d t + V aa d 2 x ) G , aQ (x, t) = S(x)S(t) (85) 

a = 1, dim K. 

We notice that the pair (a, x) is simply a label of the <fi 
field. We can choose another label (a, x) to eliminate the 



velocities and further simplify the above equation. The 
two sets of labels are related by 



(a,x) = {a,V aa x) 
Now the equations for Go. aa become 



(86) 



+ di)G , aa (x,t) = S{x)6(t) (87) 



a = 1 , . . . , dim K 

The correlation function in the presence of the boundary, 
G a b{x\,X2,t) = (^4> a (xi, t)(/>b(x2, 0)^>, satisfies according 
to © 

B T G(0,x 2 ,t) = G{x!,0,t)B = 0, 

and 

(-K ac d Xl d t + V ac d 2 Xl ) G cb {x 1 ,X2,t) 



8{xi - x 2 )5(t) 



(89) 



for xi,x 2 > 0. In the diagonal basis and in terms of the 
new label (a,i), the above equation becomes 

^- (-tsd^dt + a|J G(xi,x 2 , t) = S( Xl - x 2 )S(t) (90) 

for x\, x 2 > 0. 

The boundary condition B T <f> = can be rewritten for 
B in the form (^3|) as 

k 

M0 + ,t) - -Y / ( T oUfo+ k {0 + ,t), a=l,...,k. (91) 

b=l 

It simply connects the right moving fields 4> a to the left 
moving fields 4> a +k, a = 1,--,^, through an orthogonal 
matrix Tq. With this understanding, we find that G a b is 
given by 



G ab (xi,x 2 ,t) = (^4> a (xi,t)4>b(x 2 ,0) ^ 
= (g ( Xi - £ 2 ,t) - G (S:i + i 2 ,t)T) (92) 

V / ab 

whereT.^ ^ 

Note that by restricting to boundary conditions de- 
scribed in terms of the real bosonic fields, we are ignoring 
possible symmetries not present at the Abelian if-matrix 
level. For instance, if there are two incoming v = 1 edges 
described by Fermi fields ipi, there could be a unitary 
U{2) rotation at x = rather than the orthogonal ro- 
tation described above. As an example of the meaning 
of the orthogonal matrix Tq, consider the case of tun- 
neling between v = 1/3 and v = 1 states discussed in 
section IV. The rescaled basis is K = diag(l, 1, —1, — 1), 
and the tunneling operator is exp(i((f>i — V3(j) 2 )) + h.c. 
The matrix Tq for strong tunneling is 
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(93) 



which is the same as the matrix mapping incoming quaai- 
particles to outgoing quasiparticles in Sandler et alEB. 
For instance, two incoming electrons on the v = 1 edge 
become one outgoing electron on the v = 1 edge and 
three charge e/3 quasiparticles on the v = 1/3 edge. 

From the symmetry of the equation for Go, we see that, 
as a function of x, the matrix function Go(i,t) satisfies 



TG (x,t)T = G (-x,t) 



(94) 



B T , we can check that the above 
Certainly, G also satisfies the equation 



Using (||) and B T T 
G satisfies 

To obtain the correlation function in the original basis 
and in terms of the original labeling (a, x), we need to 
start with the explicit form of Go: Go(x, t) — — hi{Y,^x — 
t). After replacing the label (a,x) by (a,x/\v a \), we find 
(no summation over repeated indices) 



G a b(xi,X2,t) = (^(f) a (xi,t)(f>b(x2,0) 

= - ln(CE 3 ) aa (xi/\v a \ - x 2 /\v b \) - t)5 ab 

+ ln((£ 3 )aa(xi/M + x 2 /\v b \) - t){T T ) ab . (95) 

If K, V and the boundary condition B in the original 
basis are given by 



K = W I Z 3 W, V = W T V iiag W, 



B-W 1 



I 
1 o 



then the correlation function in the original basis can 
be obtained from the transformation <f> — > W^ 1 ^, K = 
S 3 -» (W T )Z 3 W, V = V diag -» (W T )V diag W, B = 



J 



W 1 



T ( T J T ] , and G -» W- X G(W 



T\-X. 



G ab (xi,x 2 ,t) = 



^(^ 1 ) QC ln((E 3 ) cc ( 



\Vr\ 



X2 

\v d \ 



t)T cd {W T r db \ (97) 



The above result for the correlations of cf> allows us to 
calculate the correlation functions of vertex operators 0\ 
by exponentiation. Consider an operator V n — e m 
Far away from the boundary (x 3> \v\t), the operator 
has a correlation which is determined by Gq only: 



where 



(V n (x,t)V n (x,Q)) ~ 1A 9 " 
g n = n T An 



(98) 
(99) 



and A = W 1 (W T ) 1 . We can write W in a form 

w =(*o 4A '* 

where Rr^l € O(k), and _B st is the boost matrix of form 



(100) 



B st = exp ( ~ T 



(101) 



In this parameterization of W , A depends only on i?^ 
(or the fc-dimensional matrix b): 

A = K- / \Bj t 2 (Kf /2 )- 1 . (102) 

Near the boundary (x -C the correlation has a 

different algebraic decay 



with 



(V a (0,t)V a (0,0)) ~ 1/i 9 " 



fl£ = n T A b „ 



where from equation ( p7|) 

A fc = IU- 1 (/-T)(IU T )- 1 . 
The above can be rewritten as 

A b = 2K- 1 B(B T AB)- 1 B T K- 1 



(103) 
(104) 
(105) 
(106) 



(gg) Since it is invariant under B — > BU, (106) is valid for all 



B, not just the form B = W T \ „ T 



C. Boundary conditions compatible with periodicity 
conditions 



In the above discussion of the boundary condition (|7fj), 
we have not considered the problem that this condition 
violates the periodicity of the fields <f>. We need to take 
into account the periodic nature of <j) field: 



2-KK~ l -a c . 



(107) 



It is clear that the boundary condition B T <j) = is not 
consistent with all the periodic conditions of the <j> field. 
Since <f> and 4> + 27r-ftT _1 n c are equivalent, if B T cf> = 
is allowed, then B T (<f) + 2TrK~ 1 n c ) = should also be 
allowed. That is, we need to generalize the boundary 
condition to at least 



B T § = 27rB T K- 1 n c , n c e f c 



(108) 



One technical way to understand what has been done 
in the previous section is that we have only considered 
boundary conditions for the neutral excitations created 
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by d x <j>- (Here "neutral" does not mean electrically neu- 
tral, but rather conserving the zero mode of the bosonic 
theory.) In addition to these neutral excitations, there 
are also charged excitations created by vertex operators 
V n = e mc where n c is a vector in the E-lattice T c . Since 
the vertex operators are the primary fields of the theory, 
it is the vertex operators which we expect to have scale- 
invariant boundary conditions, rather than th e bo sonic 
fields 4>. The generalized boundary condition (112) can 
also be written as 



i, ner. 



where the rows of B are basis vectors of T b , or 



T B = Latt(S). 



(109) 



(110) 



Strictly speaking, it is the boundary condition of the 
normal- ordered exponential which is conformably invari- 
ant, and the normal-ordered version of ( 109 ) has oo 
rather than 1 on the right-hand side. 

To gain a better und erstanding of the generalized 
boundary condition ( ]10§| ), let us consider a physical re- 
alization of the termination of the edge. We start with 
an edge described by K on (—00,00). We then add the 
following potential term on (—00, 0): 



ner f 



C n cos(n 



(111) 



where the /c-dimensional lattice Tb is a sublattice of r c 
The vectors in Tb satisfy 



n^ii" -1 !!' = 0, n,n'er, 



(112) 



and C n > are very large, so that the potential con- 
sistently pins <f) to the potential minima, and opens an 
energy gap in the region (—00, 0). Such a potential leads 
to the boundary condition 



B T (j) = 2nn, n S Latt(7fe X fe) 



(113) 



From the above discussions, we can draw two conclu- 
sions. First not all B matrices are consistent with the 
periodicity properties of cf>. To specify a valid termina- 
tion of an edge, B not only must satisfy (|72|), the rows 
of B must also be in the r c lattice, or 



Latt(.B) c f c 



(114) 



Second, two B matrices, B\ and B2, give rise to the same 
generalized boundary condition if 

B x = B 2 M, M e GL(k, Z) (115) 

Such a pair of B matrices are regarded as equivalent: 

B x S B 2 (116) 



Note that the above equivalence relation for generalized 
boundary condition is stronger (i.e., has smaller equiva- 
lence classes) than the equivalence relation B\ ~ B 2 de- 
fined in ( [73| ) for the simple boundary cond ition B T <f> = 0. 
The equivalence classes ( defined by jllo\ ) ) of the B ma- 
trices that satisfy ( J7l[ ) and (H4) label different termina- 
tions (or fixed points) of the edge. 

Now the question is what are the allowed vertex oper- 
ators on the boundary. A boundary vertex operator has 
the form VJ = e d where 1 is vector in a A: dimensional 
lattice T qB (called the boundary quasiparticle lattice). 
To determine T qB we note that Vj changes one boundary 
condition B T 4> = 2nn to another B T = 2n(n+B T K^l). 
Thus in order for B T 4> = 27r(n + B T K^X) to be an al- 
lowed boundary condition, B T K _1 l must an integral vec- 
tor. Also, we require that only shift the combination 
B T <f>. In particular, does not shift the combination 

\ B . This le 
for any 1, 1' in T 
to determine V 



r T 0, 1' S T qB . This leads to the condition fK' 1 !' 



qB '■ 



H' = 

The above two conditions allows us 



r — 

1 qB ~ 



(117) 



L&tt{KB{B T By 1 - ^B(B T By 1 B T KB(B T B)- 1 ) 



The scaling dimension of is given from (106) by 



h b (l) = 1 K B(B AB) B K I 



(118) 



In the above discussion of termination points, we have 
ignored any symmetry properties and the related selec- 
tion rules. In particular, the boundary condition charac- 
terized by B may not conserve the electric charge. As a 
result, a boundary vertex operator may not carry a defi- 
nite electric charge. In order for the termination labeled 
by B to conserve the electric charge, we must require the 
B matrix to satisfy 



B T K- 1 q = 0. 



(119) 



In this case, we find that B T <f>, the fields that are about 
to be set to a constant, commute with the electric charge 
density operator p e . For charge-conserving termination 
points, the electric charge of a boundary vertex operator 

A" 



V, b is found to be 



(120) 



The condition ( Jl 19| ) ensures the vertex operators of form 
e'^lnserB are all neutral, so that they can be set to one 
without violating the charge conservation. For a general 
B, the above vertex operators carry nonzero charges and 
setting them to one violates the charge conservation. 

For a general boundary condition B, the charge is not 
conserved, but some other quantities may be conserved. 
On the edge there are 2k conserved currents (at low en- 
ergies) j a — d t (j)a, one for each branch. Near the bound- 
ary k combinations of the 2k conserved currents remain 
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conserved. These k combinations are given by B T dt4>- 
Thus a boundary operator V b l carries k definite com- 
bined charges: 



Q = B T K-H. 



(121) 



To determine the stability of a fixed point, we also 
need to know which boundary operator VJ 6 can appear 
in the boundary Hamiltonian. First let us discuss the 
corresponding issue along the edge. Along the edge, the 
lattice r c label all the mutually local operators. Some 
carry fcrmionic statistics, and thus are not allowed in 
the edge Hamiltonian. Only the subset described by T c q 
can appear in the edge Hamiltonian. (r c o is formed by all 
the bosonic operators in T c .) If the charge is conserved, 
we further require the operators in r c o to be neutral. 

On the boundary, only a subset of the boundary op- 
erators can appear in the boundary Hamiltonian. Since 
there is no statistics within the + 1 dimensional bound- 
ary, we only need to check the conservation of the k com- 
bined charges. The values of the combined charge Q ai 
a = 1, .., k allows us to determine which boundary oper- 
ators can appear in the boundary Hamiltonian, as in the 
examples of the previous section. 



D. Image-charge picture and nonchiral fields 

This section shows how the correlation functions can 
be calculated from a simple image-charge picture when 
the chiral bosonic fields are unified into nonchiral bosons, 
as is important for a number of applications. In particu- 
lar, we find the falloff of the expectation value away from 
the boundary of a vertex operator e** pinned to 1 at the 
boundary, and how the two-body correlations are affected 
by the boundary. The correlation functions of vertex op- 
erators found in subsection B are essentially quite simple: 
any correlation function of a vertex operator can be writ- 
ten as a product of exponentials of correlation functions 
of free chiral bosons. One subtlety is that after rescal- 
ing there may be more terms in these correlation func- 
tions than experimental points in the original problem, 
since fields at the same physical point become different 
points in the rescaled coordinates. Some additional struc- 
ture appears in the correlations when the chiral fields are 
combined into nonchiral fields on the half-line, as in the 
boundary sine-Gordon model. 

In practice it is useful to combine the chiral fields <j>i 
on the whole line into nonchiral bosons $i defined on the 
half- line x < 0, in cases where the fields for x > are 
the same as those for x < 0. As an example, the appli- 
cability- of the integrable boundary sine-Gordon model 
used incl to determine tunneling behavior depends on the 
mapping to the half-line. The system on the half-line can 
be understood as a classical system on the x > half of 
the (x, t) plane, so that the physics known about such 
statistical-mechanical systems with boundaries is appli- 
cable. The technical motivation is that the theory on the 



half-line will be invariant at the fixed points under all 
the conformal generators which preserve the line x = 0. 
In what follows we show that the correlation functions of 
the nonchiral fields can be understood from an "image 
charge" picture (similar to electrostatics), and that the 
tunneling fixed points can be understood as "ordinary" 
and "extraordinary" transitions on the half-plane. 

For each chiral boson field on the whole line <pj, de- 
fine the nonchiral field &j on the half-plane x < by 
$>j(z) = 4>j(z) + </>j(z), where z = t + ix. (We use imag- 
inary time t so that conformal invariance is manifest.) 
Note that if z has x > 0, z has x < and that <&j has 
both left-moving and right-moving parts. The vertex op- 
erators exp(ia$ : j) will have different behavior depending 
on whether (j) changes sign at x = 0. First, with r\j = ±1 
the sign gained by <j)j across x = 0, for nonzero a 



(e 



lQ*j(0,t)\ _ 



> = { 







V = l 

'? = -! 



(122) 



Here and in the sequel we use the normal-ordered ex- 
ponential, which has maximum value oo rather than 1. 
Also, below we will consider the case where $j(0,i) is 
not pinned simply to but to some set of values. The 
profile of the order parameter near the boundary can be 
calculated simply: 



_( 7] = 1 

~ \ (2x)-" 2 T] = -1 ' 



(123) 



is the simplest case of the image-charge idea 
' 28 v a correlation function of n nonchiral fields 



The a 
of Card; 

on the half-plane is expressed as a correlation of 2n chi- 
ral fields on the full plane. For the boundary conditions 
we are considering, the full-plane correlation functions 
are known from (K7n, so the half-plane correlation func- 
tions of $j can be determined. The two-body function 
shows different scaling along the boundary from that in 
the bulk: (here index j suppressed) 

^ e ia*(zi) e -i^*(2 2 )^ 

_ / e M<«*i)+<K«i)) e -iW(*0+<«»2))\ 



5(a-/3)( T 



4 X!X 2 

\zi-z 2 \^z 1 -z 2 \* 

\Z!~Z 2 \ 2 ^ 

(2x 1 )o(2x 2 )P\z 1 -z 2 \2<*P 



rj = 1 
'? = -!• 



(124) 



For example, in the rj — 1 case, the equal- a; correlation 
falls off as (ti — t2)~ Aa for t% — £2 x, while far from the 
boundary (t\ — t2 < ^.x) the falloff is only as (ti — t2)~ 2a , 
i.e., with the bulk scaling dimension. For r\ = — 1 the 
correlation along the boundary is constant at long dis- 
tances, with leading correction (ti — £2)^" ■ The critical 
theory with 77 = — 1 corresponds to the "extraordinary" 
transition in statistical mechanics, where the boundary is 
ordered (the order parameter exp(i$) has nonzero expec- 
tation value) but the bulk is not, while the 77 = 1 theory 
corresponds to the "ordinary" transition. 



1G 



Some aspects of the above picture change when the 
field $(0,t) is pinned to more than one value, e.g., to 
the minima <E> = 27rrn, n £ Z of cos($/r). Now there 
is an additional average over $o — 0, ±27rr, ±47rr, ... in 
the correlation functions. The two-body correlation is 
unchanged, but the one-body correlation for 77 = — 1 is 
now 

( e «**j(z,*)) _ J (2x)~ a a — n/r,n£Z ^ (125) 
1 otherwise 

which is natural as only those operators invariant under 
the symmetry transformation <f> — > <f> + 2irr can have 
nonzero expectation values. 

VI. EDGES WITH RANDOM HOPPING 

In the %LL theory, the edge of a bulk QH state with n 
condensates is described by two symmetric n x n matri- 
ces, K and V. The integer matrix K is determined by the 
bulk QH state and is the same for all samples of a given 
edge. The positive matrix V contains non-universal ve- 
locities and interaction strengths which are expected to 
vary from sample to sample. In this section, we study 
the RG flow of the V matrix in the-presence of impu- 
rity scattering toward fixed pointsErEII which describe an 
equilibrated edge.. 

The RG calculation is given in some detail in an ap- 
pendix because there are several new features not present 
in similar treatments of the 2D classical XY modeO'EHl 
and 2p_melting£il as well as ID disordered quantum elec- 
tronso Calculations on ID quantum disordered systems 
differ from those on classical 2D systems in that quenched 
disorder is random in space but constant in time, so the 
two spacetime dimensions enter asymmetrically. The chi- 
rality of the %LL is responsible for the differences between 
our results and previous results on disordered electrons 
in ID: correlation functions in a xLL depend on x + ivt 
rather than just the magnitude x 2 + v 2 t 2 , and the oper- 
ators of interest can have nonzero "conformal spin" (dif- 
ference of right- and left-moving dimensions). One of 
the resulting RG equations disagrees with a result-pre- 
viously obtained by Kane, Fisher, and Polchinski.EI We 
outline our results before proceeding to the calculation 

In a maximally chiral edge, such as IQHE edges or 
v = 2/5, whether a given impurity operator (i.e., type of 
impurity scattering) is relevant depends only on K, not 
on V. For IQHE edges and also for the main-sequence 
chiral FQHE edges v = 2/5, u = 3/7,..., there are 
relevant impurity operators which decouple the charge 
mode from the neutral mode(s). The charge mode must 
decouple and the neutral mode velocities must equili- 
brate for the system to flow to the U(l) x SU(n) fixed 
point (n = dim _fQ —where the impurity scattering can 
be "gauged away."Lrt!J The UW) x SU (n) symmetry pos- 
sessed by K for these edgesE3 is generically broken by 
V, but restored if V flows to a decoupled charge mode 



(the U(l)) and n — 1 neutral modes with identical veloci- 
ties (the SU(n)). To our knowledge it has not previously 
been shown that general initial conditions flow toward 
this fixed point for chiral edges. The charge mode ve- 
locity is not required to equal the neutral mode velocity, 
since the disorder drives V to be diagonal in a basis where 
no disorder operator couples charged and neutral modes. 
There remain operators which couple the neutral modes 
to each other; thus although to leading order in the dis- 
order strength the neutral mode velocities do not flow 
together, it seems clear that the eventual strong-disorder 
fixed point will have equal neutral mode velocities but a 
possibly different charge mode velocity. We find that the 
fixed point is only stable if the charge mode has greater 
velocity than the neutral modes. 

The main-sequence nonchiral FQHE edges v = 
2/3,3/5, . . . have similar U(l) x SU(n) fixed points. Im- 
purity scattering now can be either relevant or irrelevant, 
depending on V , and if it is irrelevant the system will not 
flow to the fixed point. For v — 2/3 KFP used a pertur- 
bative calculation for weak-disorder to find the basin of 
attraction of the fixed pointQ; this calculation is similar to 
ours, although we find a slight disagreement (Appendix 
A). The flow to the fixed point has a much more pro- 
nounced effect on some observable quantities than in the 
chiral case: away from the fixed point, conductance and 
tunneling properties are nonuniversal. The differences 
between chiral and nonchiral edges result because scal- 
ing dimensions of vertex operators are independent of V 
(fixed by K) in chiral edges but depend on V in nonchiral 
edges. 

The edge theory of each daughter state of v = 1 
in the hierarchy is essentially the same as that of the 
corresponding daughter state of v = 1/3,1/5,.... Ev- 
ery principal hierarchy state, chiral or not, with neutral 
modes parallel to each other has a single solvable charge- 
decoupled fixed point in the presence of disorder. Edges 
with neutral modes traveling in both directions, such as 
v = 5/3 and v = 5/7, can havjg-,infinitely many fixed 
points of several different types. til The RG shows how 
for all of these fixed points the charge mode decouples, 
while the neutral modes can reach different equilibria, 
with consequences for tunneling experiments. The fixed 
points not solvable by the KFP method have disorder 
operators which frustrate each other at the fixed point. 

The xLL action in imaginary time for a clean edge 
of a QH state characterized n by the matrix K contains 
n = dim K bosonic fields cf>i'B 

S° = ^~ J dxdt [ K ij9x(t>idt4>j + V i; jd x (f>id x (j)j], (126) 

where the sum over repeated indices is assumed. K is 
a symmetric integer matrix and V & symmetric positive 
matrix. K gives the topological properties of the edge: 
the types of quasiparticles and their relative statistics. V, 
the velocity matrix, is positive definite so that the Hamil- 
tonian is bounded below. The electromagnetic charges of 
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quasiparticles are specified by an integer vector t and the 
filling factor is v = ti(K~ 1 )ijtj. 

Now a term representing quenched random impurity 
scattering is added to the action: 



Si 



dxdt [£(a:)e im ^ + C(x)e~ m ^^ 



(127) 



Here £ is a complex random variable and ({£(x)£* {%'))) = 
£><5(ir — x'), with D the (real) disorder strength. The 
integer vector m describes how many of each type of 
quasiparticle are annihilated or created by the operator 
O m = exp(imj(j)j). For a real system all charge-neutral 
scattering operators rrij are expected to appear, but most 
of these will be irrelevant in the RG sense. The condition 
for charge-neutrality is ti(K~ 1 )ij r mj = 0. The random 
variables £ m for different scattering operators O m may 
be uncorrelated or correlated depending on the nature of 
the physical impuriti es ca using the scattering. 

The clean action (126) is quadratic and hence does 
not fl ow under RG transform ation s. Adding impurities 
(|127|) causes the V matrix in (126) to flow, and in some 



cases .the, flow is to a new type of strong-disorder fixed 
point ntl Here we will quickly review the diagonalization 
of the clean action to find the eigenmodes and their 
velocities v i: and then find the RG flows for two- mode 
edges with a single impurity operator. Then a general 
edge with several modes and impurity operators is con- 
sidered. 

Let M\ be some matrix which brings K to the pseudo- 



identity I n 



,-: K- 



MiI n+ , n _Mi . Then V can be 



brought to a diagonal matrix V D = M 2 T Mi T VMiM 2 , 
where M 2 is an element of the group SO(n + ,n~) so that 
M 2 Mi still takes K to the pseudo-identity. The point of 
these transformations is that the action is now diagonal 
in the basis <fi — (MiAfa) -1 ^, so the correlation functions 
are simple: 



( e i<Pj (x,t) e -i^ (0,0) \ _ e (<t>j {x x t)<t>j (0,0)) - (<t>j {0,Q)<t>i (0,0)) 



oc (x ± ivjt) 1 



(128) 



where the sign depends on whether cf>j appears with —1 
or +1 in I(n + ,n~). The vertex operator O m described 
by the integer vector m has correlation function 



l e imj<jij(x.t) e 



-im.j(f>j(0,0)\ 



= ]]_(x±iV j ty 



(129) 



with mjifij = Cj<j)j. The total scaling dimension of O m is 
A(m) = c |/2j which is bounded below by K(m)/2 = 
mK~ 1 m/2. The impurity term Si containing O m with 
a random coefficient is relevant if A(m) < 3/2; the cor- 
responding marginal value for a uniform coefficient is 2, 
and for a S- function coefficient 1. 

Appendix A calculates the change in the correlation 
function of O m under an infinitesimal RG transforma- 
tion induced by the impurity term Si. Here we find the 



change in the underlying V matrix required to produce 
the new correlation function. The K matrix is unchanged 
as it is "topological" (it does not enter the Hamiltonian). 
The V matrix flow has a simple interpretation, valid for 
any number of edge modes traveling in cither direction. 
Each impurity operator m drives V to become diagonal 
in the basis with m an eigenvector. This automatically 
minimizes the scaling dimension of O m in a nonchiral 
edge. In cases where there are more impurity operators 
than independent eigenvectors, so that not all impurity 
operators can simulatcnously be eigenvectors, the impu- 
rity operators frustrate each other. 

Both v = 2 and v = 2/5 have a single K(m) = 2 
operator which is always relevant. In the basis ei = t, 
e 2 = m, 



K 1 ' 2 



1/y/v 
1/2 



V = K 1/2 R 



Vi 







V 2 



(130) 



with R a two-by-two rotation matrix by some angle 9. 
Note that all V are obtained by considering 9 in the in- 
terval [Q, 71") . Now a = 2sin 2 and (3 = 2 cos 2 (9 are the 
exponents appearing in the correlation function of the im- 
purity operator: (O m (x, t)Oj n (0, 0)) — (x + ivit)~ a (x + 
iv2t)~P, a + /3 = 2. Then to first order in disorder 
strength, the diagonal velocities v\ and v 2 are unchanged, 
and (Appendix A) 



da 

m 



8nDa(3 



(vi - v 2 )vi a 1 v 2 P 



Since da = 2sin(29)d9 = 2y^pd9, 

d9 _ 4irDsm(29)viv 2 
dl ~ (vi-v 2 )vi 2sin20 v 2 2cos2 < 



(131) 



(132) 



There are two fixed points of this equation, with 9 = 
stable and 9 = tt/2 unstable for m > v 2 , and vice versa 
for vi < v 2 . The stable fixed point always corresponds to 
neutral mode velocity less than charge velocity (Fig. ^i) . 

We can summarize the effect of the disorder opera- 
tor in the comoving case simply: it rotates V so that 
m becomes an eigenvector. Since m is neutral (mii' _1 t 
= 0) the other eigenvector is driven to the charge vec- 
tor. The idea that impurity operators drive V to make 
themselves eigenvectors is quite general. The case of two 
countermoving modes (e.g., v = 2/3) with a K(m) = —2 
disorder operator is similar in form. The rotation matrix 
R in (130) is replaced by a boost matrix B, 



B 



( coshr sinhr 
I sinh r cosh r 



(133) 



and the exponents in the correlator are a = 2 sinh r, 



(3 = 2 coslr r, a- -- 
equation for r is then 



-2, a + /3 = 2A(m). The flow 
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dr 4ttD sinh(2r)v + v_ 



(134) 



Here u+ and u_ are the (positive) velocities of the right - 
and left-moving modes. Now there is only one fixed point, 
at r = OlFig. 0b), which is the solvable fixed point found 
by KFP.Q 



(a) 




(b) 



X = — oo 



-> • <- 

x = 



X = oo 



FIG. 6. Schematic RG flows for two-mode (a) chiral and 
(b) nonchiral edges. The same idea applies to main-sequence 
edges with more than two modes: the chiral case has one 
stable (« c > v n ) and one unstable (v c < v n ) fixed point, while 
the nonchiral case has one stable fixed point independent of 



Now we consider a general edge with several modes and 
impurity operators. To first order in disorder strength, 
the effects of each impurity operator add independently. 
Scaling dimensions A(m) of vertex operators O m are in- 
dependent of V in chiral edges, depending only on K: 
2A(m) = mK = K(m). In a nonchiral edge this 
holds as an inequality: 2A(m) > K(m), with equality 
only if V is diagonal in a basis with m an eigenvector. 
Since most experimentally relevant quantities are deter- 
mined by scaling dimensions, it is useful to isolate which 
parts of V affect scaling dimensions. The matrix M 2 
used above to diagonalize V while preserving the pseudo- 
identity /„+ „- is an element of SO(n + ,n~). It can be 
decomposed M 2 = BR into a product of a 'boost', a sym- 
metric matrix, and a rotation, an orthogonal matrix, each 
elements of SO(n + ,n~ -)B The n(n+ l)/2 free parame- 
ters in the symmetric positive matrix V are now taken 
as n eigenvelocities (the elements of Vd), n + n~ "boost 
parameters" (which correspond to interactions between 
oppositely directed modes), n + (n + — l)/2 rotation pa- 
rameters between right-movers, and n~(n~ — l)/2 ro- 
tation parameters between left-movers. Only the n + n~ 
boost parameters affect scaling dimensions because, in- 



troducing the matrix Ay via 2A(m) = HijAjjinj, 

A = M 1 BRR T B T M 1 T = M 1 B 2 M{ r . (135) 

For each pair of comoving modes appearing in the cor- 
relation function, there is an infinitesimal change in the 
rotational part R, and for each pair of countermoving 
modes, there is an infinitesimal change in the boost part 
B. 

Given an impurity operator O m and initial V{£), we 
need to find V(£ + d£) which gives the changes to the 
correlation function calculated in Appendix A. However, 
there is an important issue not present in the two-mode 
case: there are more free parameters in V than expo- 
nents in the correlation function, so V is not uniquely 
determined without additional assumptions. We assume 
that each term in Appendix A coupling two modes affects 
only the components of V between those two modes. We 
will also write the flow equations for the components of 
V directly, rather than introducing a parametrization as 
we did above in terms of 6 or r, because for multi-mode 
edges such parametrizations become quite complicated. 

Suppose that under an infinitesimal change d£ the V 
matrix changes from (MiM 2 )~ 1T V D {M 1 M2)~ 1 to 



V(£ + d£) = {M l M 2 N)- 1T V D {M 1 M 2 N) 



(136) 



Note that V does not flow if the eigenmodes have the 
same velocity, so that Vd is a multiple of the identity. 
Here N is an element of SO(n+,n—) differing only by 
order d£ from the identity. The fields which diagonal- 
ize V(£ + d£) are <j>' = [M^N)' 1 ^ = N^ 1 ^. Now let 
Cj be the components of the disorder operator in the <f> 
basis which diagonalizes V(£): rrii4>i — Cj4>j. The expo- 
nents appearing in the correlation function, whose flow 
is calculated in Appendix A, are a = ci 2 , (3 — c 2 2 , etc. 

Now we are ready to construct N. For each pair of 
countermoving modes i and j, N has an infinitesimal ro- 
tation angle dOij = —dOji, and for each pair of comoving 
modes i and j, a boost drij = drji. For an edge with two 
right-movers and one left-mover, e.g., 



1 


-dd 12 


dr 13 




1 


dr 23 


dr 13 


dr 23 





N 



Then for a disorder operator O r 
strength D, 



(137) 



with projections Ci and 



Airci Cj Vi Vj D d£ 

[Vi- VjlUi ^ 

AirciCjViVjD d£ 



(138) 



Note that the Cj can have either sign and must change 
sign in the vicinity of a fixed p oint , as above in the two- 
mode case. Substituting N in (136) gives the desired RG 
flow equation for the components of V. The condition 
for a fixed point is simple: N t VdN — Vd- Hence for 



19 



any chiral edge (N a pure rotation), if all velocities are 
the same the system is at a fixed point. 

The simplest multicomponent edges are the IQHE 
edges v — n and chiral main-sequence edges v = n/(2n + 
1). The behavior of these is similar to the v = 2 case dis- 
cussed in detail above: the stable fixed point has charge 
and neutral mode velocities v c > v n . The neutral mode 
velocities are expected to equalize while the charge mode 
remains different, because there are always hopping oper- 
ators connecting the different neutral modes, while there 
is no hopping connecting the charge mode to the neutral 
modes. As a result there is no stable fixed point unless 
the neutral velocities are the same. 

For the v = 3 case we can demonstrate the similarity 
to the v = 2 case by an explicit calculation. Parametrize 
the rotation part R of V as R = R23 {8i)R\ 2 {9 2 ) -^23(6*3), 
where Rij (0) is the rotation by 9 in the i — j plane; then 
in a basis with ei = t, the charge mode is decoupled if 
9 2 = 0. The flow equation for 2 , with v 2 = v 3 = v n , is 

d6 2 —4:nv c v n sin(2#2) 
d£ (v c - v n ) 

cos 2 fli(4£>i + D 2 + D 3 ) 

/^y 2 cos 2 6\ sin 2 &2y 2(sin 2 #i+cos 2 0± cos 2 #2) 

3 sin 2 0i(£> 2 + £> 3 ) 

4 Uc 2sin 2 81 sin 2 S 2Uji 2(cos 2 Sx+sin 2 8 ± cos 2 9 2 ) 



(139) 



where -D1...3 are the strengths of th e 3 hopping operators. 
This is of the same form as (132) since the quantity in 
brackets is clearly positive. Henceforth we will not write 
out the flow equations but just discuss the qualitative 
behavior. 

The nonchiral main-sequence edges v = 2/3, 3/5,... 
have solvable stable U(l) x SU(n) fixed points if all the 
neutral mode velocities are equal and the charge mode is 
decoupled. The charge mode does indeed decouple in the 
perturbative RG equations since each impurity operator 
reduces the scaling dimension of t toward its minimum 
v, which is only attained when the charge mode is de- 
coupled. Once the charge mode is decoupled, the neutral 
modes behave exactly as in the chiral case, except that 
the fixed point is stable even if the neutral mode velocity 
is greater than the charge mode velocity. 

In the edges v — 5/3 or v — 5/7, which have neutral 
modes in both directions, there are an-jpfinite number 
of possibly relevant impurity operators.Eil However, near 
each of the possible fixed points there are only three rel- 
evant operators. Unlike the case with all neutral modes 
in the same direction, to first order in disorder strength 
there is no stable fixed point for V, even if the neutral 
mode velocities are equal. The mathematical difference 
is that, now that the modes move in opposite directions, 
the impurity operators cause infinitesimal boosts which 
do not disappear when the velocities are equal, unlike 
infinitesimal rotations. The fixed points in these edges 
have marginal operators not present for the KFP-type 



fixed points and seem to be of a different type, and their 
strong-disorder behavior and stability is not well under- 
stood. 

There are several four-component edges, such as v = 
12/17 and v — 12/31, which have solvable, SU(n) sym- 
metric fixed points of the KFP type as well as fixed points 
with marginal operators similar to the v = 5/3 and 
v = 5/7 edges. The picture of equilibration differs de- 
pending on the fixed point. The solvable SU(2) x SU{2) 
fixed points found in these edges can have all three neu- 
tral mode velocities different (since there is a basis where 
no relevant hopping operator couples neutral modes), 
while the SU (3) fixed points have to have two neutral 
mode velocities equal in order to be solvable. However, 
as discussed in the introduction the fis^d point accessi- 
ble by the composite-fermion approaclJIS is the one with 
marginal operators present, even though this is the onijj 
one of the three fixed points whose stability is doubtful. EJ 



VII. CONCLUSIONS 

The main result of this paper is the analysis for 
a general Abelian edge of a large class of tunnel- 
ing/termination fixed points. We find the conditions on 
allowable boundary conditions and the resulting corre- 
lation functions for vertex operators at the fixed points. 
We find that the requirement of unitary time evolution 
imposes conditions on when a quantum Hall edge can 
terminate at a point or join smoothly to another edge. 
Two edges can join smoothly only if they have the same 
number of modes, and there are additional statistical 
restrictions: for example, v — 1 can join smoothly to 
v = 1/9 (provided a source of charge is present), but not 
to v = 1/3. We find the lattice of allowed boundary op- 
erators, both charged and neutral, for a given boundary 
condition. 

The correlation functions in the presence of a bound- 
ary can be understood through an "image-charge" pic- 
ture, and take a simple form which becomes more com- 
plex when the original chiral bosons are combined into 
nonchiral fields as for the boundary sine-Gordon model. 
The two-point correlation function of a vertex operator 
can show different behavior (a change in scaling dimen- 
sion) as the two locations are moved toward the bound- 
ary. 

We solve a simplified model of how interactions 
between electrons on different edges affect tunneling 
through a point junction. Most current tunneling exper- 
iments involve multiple contacts between the two edges 
and variable interaction strengths along the edge, but the 
basic result that interactions can complicate the identifi- 
cation of the effective tunneling charge should still apply. 
The model we consider maps onto the solvable model of 
Fendley et al. with a continuous effective filling fraction, 
giving a physical realization of this result beyond the dis- 
crete cases v e g = 1,1/2,1/3 known previously. 
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We have used a perturbative RG approach to find the 
flows of random hopping operators in a general quantum 
Hall edge. Two broad conclusions from this treatment 
are that with no initial assumptions favoring one type 
of state over another, significant differences emerge be- 
tween non-principal states and principal states, and be- 
tween states with all neutral modes in one direction and 
states with neutral modes in both directions. In states 
with neutral modes in both directions, different hopping 
operators compete to drive the system to different fixed 
points. This may explain why hierarchy quantum Hall 
states with neutral modes in both directions are much 
more difficult to observe than those at the same level of 
the hierarchy with all neutral modes in one direction, and 
why non-principal states are also rarely seen. We find 
that charge-neutral separation is a general feature when 
disorder is relevant, in both chiral and nonchiral edges. 
A specific result for the v = n IQHE states and chiral 
main-sequence FQHE states is that the charge-neutral 
separated fixed point is only perturbatively stable when 
the charge mode is faster than the neutral mode, as for 
strong Coulomb interactions. 

X.GW acknowledges support by NSF under the 
MRSEL Program DMR 98-08491 and by NSF grant No. 
DMR 97-14198. 



APPENDIX A: PERTURBATIVE 
RENORMALIZATION GROUP 

This appendix uses the perturbative RG technique to 
study the effect of impurity scattering operators on the 
velocity matrix V in the chiral-Luttinger-liquid action 
(|J) . The K matrix does not flow and thus remains an in- 
teger matrix, which follows directly from the fact that K 
does not enter the Hamiltonian (i.e., it is "purely topo- 
logical" ) . The constancy of K to leading order in D will 
be explicit in the results obtained below. In order to cal- 
culate the changes in V under a change in the cutoff, we 
expand a correlation function to first order in the disor- 
der strength, then show that the terms proportional to 
the disorder strength can be interpreted as infinitesimal 
changes in the matrix V. 

The real-space calculation is similar .-to previous RG 
calculations on 2D classical modelsEjEil and ID elec- 
tronsEl The differences arise from the chiral nature 
of quantum Hall edge states. As an example, consider 
the correlation function of a vertex operator O m — 
exp(im,j(pj) in a nonchiral edge with one mode in each 
direction: 

(Omir^OKrz)) <x(x + iv + t)~ a (x - iv-t)~^ (Al) 

with a — [3 — K(m) an even integer, r, = (xi,ti), x = 
X2 — X\ and t = t% — t\, and (v+,v—) the velocities of 
the right and left moving modes. Unless K(m) = and 
v + = u_ , the correlation function has a phase as well as 
a magnitude. 



First we treat the case of an edge with two modes, 
either parallel or antiparallcl, and then show how the 
flows for an edge with more than two modes follow with 
no further computation. 

The correlation function of an operator O n , expanded 
to first order in the disorder strength, is 

/giTty^riJg-ira^fo)^ _ Uirij^in) e -inj<j>j(r 2 )\ x 
(1 - fdr 3 dr 4 la(x 3 )e(x4)(e m ^^e- m i^) }) 



(A2) 



+ / dr 3 dr 4 [£,(x 3 )C (x 4 )x 

Jiij (n) „-inj 4>j (r 2 ) imj 4>j (r 3 ) -imj 4>j (r 4 ) \ 



Now carry out the disorder average {£(x)£* (%')} = 
D8{x — x') and consider the term with four correlation 
functions. Introduce R = (r3 + r^)/2, r = r 4 — r3. Only 
configurations where the internal points r3 and r 4 are 
near each other (i_e. , separated by the cutoff a) contribute 
to the RG flows E3 At this point assume for convenience 
that we are calculating the correlation function of the 
disorder operator itself (m = n). The symbol P12 de- 
notes (x + iv + t)~ a (x — ii)—t)~°, x = X2 — x\, t = t% — t\. 
Because a — (3 = K(m) is even, P12 = f^i- The last term 
in (|A2Q is now 



D / dXdTdt[ 



-P12-P34-P14-P23 J 
-P13-P24 



(A3) 



The integrand is P12-P34 exp(5i 3 + g 24 - 9u ~ 923) ~ 
-Pi2-P34exp(r • Vr(s(i-i - R) - g(r 2 - R))). Here g a b = 
— \og(Pab) = alog(x + iv + t) + (3log(x — iv-t). Disor- 
der fixes x = in r = (x, t) so the exponential is 
exp(t&r(s(n-R)-fl(r a -R))) ~ 1 + t 2 {d T g(ri - R) - 
drg{ r 2 — R)) 2 - Hence we need to evaluate the following 
integral: 



D J dX dT dt P 12 P 3A [l + t 2 { 



iav + 



X — x\ + iv+(T — t\) 



X — x\ — iv^(T — ti) X — X2 + iv+{T — £2) 

+ i_P^ s2l 

X - x 2 — iv-(T — t 2 ) 



(A4) 



The constant term in the integrand cancels the leading 
term in the partition function in (|A2Q. When the square 
is expanded, products of denominators at the same point 
will cancel infinities arising in the calculation of products 
of denominators at different points, leaving a finite an- 
swer. The change of cutoff in the t integral will yield the 
RG equations at the end. 

First consider the integrals of the a 2 ,/? 2 terms. After 
rescaling the time variables by v + , the a 2 integral is 



dXdT\ 



1 



1 



/ 



dXdT\ 



X + iT — x\ — it 1 X + iT — X2 — iti 
1 , 



(A5) 
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with z = X + iT,2w — x 2 — x\ + it 2 — it\. But this 
integral is not uniformly convergent at 00 and thus not 
well-defined: for example, if to = i the integral is 2tt 
if the X integration is done first, if the T integration 
is done first, and 7r if the integration is done in radial 
coordinates. We believe that the appropriate value of 
the integral (A.5) is 0, because in Minkowski space (real 
rather than imaginary time) the corresponding integral 
has integrand (x + t — x\ — ti)~ 1 (x + t — X2 — fe) -1 and 
is unambiguously zero. Also, zero is the only value con- 
sistent with the fact that the randomness can be rotated 
away at the KFP fixed point, since at that point the RG 
flow of the velocity should be independent of the disorder 
strength. 

The a(3 terms, give the renormalization of the scaling 
dimension A, are proportional to 



I 



dXdT\ 



1 



1 



c_d + c+d- 
c± = X — x\ ± iv±(T — <i), 
d± = x 2 - X ± iv±(t 2 - T). 



(A6) 



First do the dX integral as a contour integral. The poles 
of the first term are at toi = X\ + ii)-(T — ti) and w 2 — 
x 2 + iv+(t 2 — T), and the integral vanishes unless the 
poles are on different sides of the real axis (likewise for 
the second term). Thus T ^ [ii, £2] and we are left with 



I = 



dT 



12 



dT 



2m 



2m 

V2 



Vi = 
V2 = 

Hence 
I = 



x 2 
x 2 



2tt 



x\ + iv + t 2 + iv-ti — i(v + + V-)T, 
x\ — iv + ti — iv-t 2 + i(v + + v-)T. 



(A7) 



v 



a 



± 



x 2 

do = X 2 

= x 2 - 



log 



aV ffli a, a 9 



(b+b-) 2 

x\ ± i(v + + v-)oo + iv + t 2 + iv-ti 
x\ ± i{v + + v-)oo — iv + t\ — iv-t 2 , 
x\ ± iv±(t 2 — t\). 



(A8) 



The finite part of the result is independent of the order 
of integration, eve n th ough (A6) is superficially even less 
well-defined than (|A5|). The i nfin ite part of the result is 
canceled by the a/3 terms in (A4) with denominators at 
the same point. The dt integral is 



dtt 2 P^A = 2 



1 1 

dt[t 2 (l—T(-r—f] 



Hv+t —iv-t 

roo 



v + a vJ j 



dtt 2 - (a+p) . (A9) 



The dependence on K(m) here is an artifact of our 
using the unfortunate convention (Al), when in fact 
the proportionality constant in front alternates sign as 
|Jf(m)| = 0, 2, 4, ... to keep the correlation function pos- 
itive when its argument is on the time axis. The effective 
scaling dimension A e ff after re-exponentiating the per- 
turbation to the correlation function is 



2A cff = 2A - 



4(47r)a/3D 



("4 



p-1 



dt 



ta+P-2 • 



(A10) 



By the usual process of changing the cutoff a — > a exp(^) 
we obtain the RG equations 



^-(3 
dl ~ (6 

dA _ 



-2A)D 
8tt£)(A 2 



K(m) 2 /4) 



/3-r 



(All) 



These equations match those found by Kane, Fisher, and 
Polchinski for the |Jf(m)| = 2 operator in the v = 2/3 
state. However, as mentioned above we find no term 
which rcnormalizes the velocities to first order in D. 

Now consider the case of two comoving modes. The 
simplest example is the IQHE state v = 2, which has a 
relevant impurity operator hopping electrons from one 
mode to the other. The correlation function P\ 2 of 
an impurity operator O n — exp(irij(f>j) has the form 
(x + ivit)~ a (x + iv 2 t)~P , with (vi,v 2 ) the velocities of the 
two eigenmodes and a + (3 — K(n) an even integer. Ex- 
panding the correlation function in the disorder strength 
and then evaluating the disorder average as before gives 
the correction to the correlation function: 



D 



dXdTdtP 12 P 34 [l + t 2 { 

if3v 2 

X — x\ + iv 2 (T - ti) 

i/3v2 ^ 

X - x 2 + iv 2 (T - t 2 ) ' 



iav\ 



X 



x\ + ivi(T - 
iav\ 



*i) 



X 

2l 



x 2 + ivi(T - t 2 ) 



(A12) 



Expanding the square gives terms with both denomina- 
tors at the same point, which cancel infinities appearing 
elsewhere in the calculation, and terms with both denom- 
inators having the same velocity which were previously 
argued to be zero (and in any event cannot cause the two 
velocities to flow toward each other, since each term only 
involves one velocity). The result is 

1 



a(3vxv 2 P 12 J dXdT[ 
' 1 



X ~ x 2 + iv 2 (T - t 2 ) 
ATta[3viv 2 Pi 2 J dtt 2 P 3 . 
vi - v 2 



X - xt + ivi{T - h) 
+ (vi ~v 2 )} { dtt 2 P 3i 



log 



x 2 - x x + iv\(t 2 - t x ) 

X 2 — 351 + W 2 (t 2 - tl) 



(A13) 



Thus a and (3 are changed but not a + (3 = K(m). The 
velocities of the eigenmodes are unaltered, and the RG 
flows are 



dD 

~dl 

da 
HI 



(3 — a 
dp _ 
~~d£ ~ 



■0)D 



(3 - K{m))D 

8nDa(3 



(vi - v 2 )vi c 



-1^/3-1- 



(A14) 



The singular denominator when v± = v 2 is acceptable 
because at v\ = v 2 , only a + (3 is well-defined, not a and 
(3 separately. 
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The generalization to a case with two modes and more 
than one impurity operator is simple: the contributions 
to the RG flow equations for the velocity matrix from 
each impurity operator add, since to leading order the 
impurity operators are independent. 

Extending the calculation to an edge with more than 
two mod es is quite simple. In the expansion of the square 
term in (A4), each pair of modes gives one term. If the 
two modes move in opposite dire ction s, the term lowers 
the total scaling dimension as in ( All ); if the two modes 
move in the same direction, the term maintains the to- 
tal scaling dimension as in (A14). Each term preserves 
K (m) separately. For explicitncss, consider the case of an 
edge with two right-movers and one left-mover, which is 
relevant to the v = 3/5 edge. Then the correlation func- 
tion has the form (x + iv\t)~ a {x + iv2t)~P(x — iv-t) 1 , 
and the RG flows for the exponents are 



8irDa"/viV- 



da 8irDaf3viV2 

d(3 8irDa(3viV2 8irD/3jV2V_ 

di (v\ — V2)V\ ol V2^V- 1 (l»2 + V-)vi a V2^V-^ 

c?7 8irDajv\V- 8irDfl'y'W2 I V- 

di (v\ + V-)v\ a V2^vJ^ (V2 + V-)v\ a V2 l3 V- 1 

Section VI shows how the V matrix flow in the x^L 
action determined by the correlation function flows found 
in this appendix has a natural interpretation in terms of 
the boost and rotation parts of the V matrix. 
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